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The dynamics of a linear system embedded in a heat bath environment and subject to non- 
Gaussian noise is studied. Higher order cumulants in coordinate space are derived and their impact 
on the dynamics and on asymptotic steady state distributions is analyzed. In the quantum regime 
non-Gaussian properties are present in the reduced density in coordinate representation which in 
energy representation exist on a transient time scale only due to symmetry. Within an exactly 
solvable model our results provide insight into mechanisms of linear detectors as sensors for non- 
Gaussian noise at high and low temperatures. 
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I. INTRODUCTION 

Non-Gaussian noise plays an important role in a va- 
riety of fields including nanoscopic solid state structures 
1 , Q , electrical circuits Q , signal and image processing 
4 1, and biological aggregates [f|. In most cases addi- 
tional thermal fluctuations are also present due to the 
interaction of the system of interest with a macroscopic 
background constituting a heat bath. According to the 
dissipation fluctuation theorem, the latter one gives also 
rise to energy exchange such that in absence of non- 
Gaussian noise the system would equilibrate during its 
time evolution. The presence of non-Gaussian random 
forces, typically stationary but out of equilibrium, drive 
the compound towards a steady state though which may 
significantly deviate from its thermal equilibrium. 

Theoretically, a systematic description for these pro- 
cesses has been developed in the context of charge trans- 
fer through mcsoscopic conductors. Within the frame- 
work of full counting statistics low frequency noise prop- 
erties for various devices have been determined 0]. For 
simple devices such as tunnel junctions, predictions have 
been verified experimentally, cf. e.g. [g, but the detec- 
tion of non-Gaussian noise signatures at finite frequen- 
cies is a challenging task. The same is true for a theoret- 
ical understanding which in addition to the noise source 
must also include the detector degrees of freedom. In the 
classical regime, progress has been achieved with Joseph- 
son junctions as threshold detectors Q. In this set-up 
the dynamics of the superconducting phase is the rele- 
vant coordinate which is subject to thermal noise and 
non-Gaussian fluctuations originating from a mesoscopic 
conductor and acting as an external driving force. Its 
random motion corresponds to that of a fictitious parti- 
cle in a tilted washboard potential such that the rate of 
escape from one of the metastable wells encodes infor- 
mation about the non-Gaussian noise An extension 
to the quantum regime including tunneling through the 
barrier is still elusive, however. 

In this paper we consider a much simpler situation, 



namely, linear dynamics in presence of non-Gaussian 
noise, a problem which allows for a full analytical treat- 
ment. However, the results are by no means trivial and 
due to the ubiquitous appearance of harmonic systems 
in all fields of science may be of general interest. In fact, 
in the context of Josephson junctions this corresponds 
to the motion of the phase deep inside a well which is 
known to provide already the correct parameter depen- 
dence of the dominating exponential of the escape rate. 
There is also another interesting facet of the subject: In 
the quantum regime the reduced density operator (ob- 
tained after integrating out thermal and non-Gaussian 
noise) is the relevant object and it is the measurement 
process which determines its appropriate representation 
in Hilbert space. We show that the presence of non- 
Gaussian properties in the steady state very sensitively 
depends on this representation such that e.g. for mea- 
surements sensitive to energy eigenstates the skewness 
(third order cumulant) is absent while it can be observed 
in coordinate space. These findings lay on the one hand 
the basis for implementing a linear detector as a sensor 
for non- Gaussian noise and on the other hand pave the 
road for the treatment of anharmonic systems. 

The paper is organized as follows: General results in 
the classical regime are presented in Sees. HU IIIII which 
will then be specified to the case of Poissonian noise in 
Sec. IIVI The linear problem can be mapped also to the 
case of free Brownian motion with Drude damping as 
illustrated in Scc.|Vl a simple model for a tunnel junction 
in series with an RC-circuit as impedance. The quantum 
mechanical problem is treated and discussed in Sees. IVII 
and lVIH Conclusions are given in Sec. IVIIll 



II. PRELIMINARIES 

A linear system with generalized coordinate q embed- 
ded in a heat bath and subject to an external force F{t) 
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can be described by the generalized Langevin equation 



q(t)+ / ds^{t-s)q{s)+LJ 2 q(t) 
Jo 



±[£(t) + F(t)]. (1) 



Here, £ is a fluctuating Gaussian force with = 

which is related to the damping kernel "/(t) via a dissipa- 
tion fluctuation theorem, i.e., (£(£)£( s )) = ~ 
s), t > s, with inverse temperature /3 = 1/k^T. 

In the sequel, we will concentrate on the simplest case 
of white Gaussian noise 7(t) = 2t<5 (t), but the gener- 
alization is straightforward (see e.g. The exter- 
nal force is taken as a random non-Gaussian variable 
F(t) = rj(t) with (rj{t)) = 0. It is assumed to be statisti- 
cally independent from the thermal force (£(t)r)(s)) = 
and (^-correlated (white non-Gaussian noise). While this 
certainly simplifies actual experimental situations, it will 
nevertheless provide insight into the major impact of non- 
Gaussian noise on the classical and quantum dynamics of 
a dissipativc oscillator. 

Accordingly, we have solutions of the form q(t) = 
Qo{t) + q^(t) + q v {t) with the homogeneous part 



qo (t) 



i* ( Cl e^ 



C2e 



(2) 



with to = \/u!q — 7 2 /4 and constants C\ , C2 determined 
by the initial conditions. The force term follows from 
the response function of the problem 

G(t = 6{t tl)e -^-t>)^Ht-t')) (3) 

Moo 



and reads 



(t) = f dsG(t-s)y(s) 
Jo 



(4) 



for the thermal noise (y = £) and the non-Gaussian driv- 
ing (y = rj), respectively. 



III. CUMULANT-GENERATING FUNCTION 

We are interested in position correlation functions of 
the oscillator beyond second order and particularly their 
steady state properties. For this purpose, it is convenient 
to work with cumulant generating functions 



Fly H 



y = t,v (5) 



Here, w(u) plays the role of a time dependent counting 
field. Correlation functions of the full solution q(t) = 
Qo(t) + + Qr](t) can then easily be calculated in a 
straightforward way. Since the noise terms affect only the 
inhomogeneity of the differential equation, the additive 
homogeneous solution go provides trivial contributions 
and will be ignored in the sequel. We thus concentrate 
on the deviations from the deterministic motion. The 
n-th order cumulant then follows from 



C qy (t n ,...,ti) = {q y {t n )...q y {t 1 )) y 

xn 

= nr 



Sw(t n ) ■ ■ ■ 8w{t\ 
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where t n > t n -\ ■ ■ ■ > t\. Cumulants for the momenta 
p y (t) = Mq y are obtained from taking the corresponding 
time derivatives. 

Now, given the linear structure of the Langevin equa- 
tion ([I} , the generating functionals for q y can be directly 
inferred from the generating functionals of the respective 
noise forces y, i.e., 



F qy [w] = F y [w] = In^/o" *»v(«)*C«)) v 
with a redefined counting field 



(7) 



w{u) = I ds G(s — u)w(s) = I ds G(s — u)w(s) . 
Jo Ju 

. (8) 

Due to the linearity (harmonic system), there is a 
linear mapping from its cumulants to the character- 
istics of the noise source via the response function. 
Accordingly, for the thermal noise one has = 
(— M7//3) J duw(u) 2 so that the only non-vanishing 
cumulant reads 



C qg (t, t') = (M7//3) / duG{t - u)G(t' -u);t>t'. 

(9) 

A compact expression is obtained for the equal-time cor- 
relation 



1 



M/3uj 2 



Alu 2 — 7 2 cos(2wt) — 27wsin(2wt) 



4M/3uj 2 uj 2 



(10) 

which asymptotically reduces to the equilibrium variance 
as expected. 

To make further progress with the non-Gaussian noise 
contribution, we specify results to Poissonian noise 

(r)(u k ) ■ ■■v( u i)) v = \a k S(ui-u 2 ) . ..6(uk-i-u k ) , k > 1. 

In an electrical circuit noise of this characteristics is ap- 
proximately generated by a voltage-driven tunnel contact 
with a voltage much larger than k^.T . Then, a corre- 
sponds to the transferred charge and A is the transport 
rate. The cumulants in position for the harmonic oscilla- 
tors follow from the known generating functional of the 
Poisson process 



F n [w}=\ j &u(c law{u) -\). (12) 

—00 

according to T q [w] = J- v [w] . 



IV. CLASSICAL CORRELATION FUNCTIONS 
IN POSITION SPACE 

With the generating functional at hand we can now an- 
alyze the lowest order cumulants in more detail. The sec- 
ond cumulant is trivial as it is identical to C q . up to a con- 
stant factor, namely, C q (t,f) = \a 2 /(2M~// /3)C q At,t'). 
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Thus, the presence of the Poissonian noise leads to an 
additional heating and it can simply be taken into ac- 
count by introducing an effective temperature T c g = 
T + Xa 2 /{2Mjk B ). 

More interesting is the three-point-correlation function 



Xa 3 



w(4w 2 - 7 2 ) 



-S-yt 



M 3 w 3 |6w 2 (2 7 2 +u; 2 ) 
2ui cos(3wt) + 7sin(3a;i) 



24w 2 



-9^5 



2 2w cos(a;i) + 37 sin(wt) 



2 7 2 



■UK 



(13) 



which yields asymptotically the stationary expression 



a 



(3) 
Qv 



lim^oo Cq v (t, t, t) with 



2Xa 3 



3M 3 w 2 (2 7 2 +w 2 ) 



(14) 



(2) 

Hence, in contrast to Cq n , the skewness becomes basi- 
cally independent of friction for weak dissipation. 

All higher order cumulants can be calculated accord- 
ingly. After exploiting (jll[) the general expression re- 
duces to 



H.TN,.,N 



AT 



N-m 



in 



m— OIIijv i—1 j — 1 



ti 



x(-l) m / due" M sin m (a;u)cos iV - m (a;u), (15) 


where in the second sum the indices Pi run over all per- 
mutations HiN of the set {1, . . . , N}. While no compact 
expression for the remaining integral is known, it can 
easily be calculated for specific values of N, m. 



The stationary cumulants C, 



(AQ 
Qt, 



limt^oo C q (*,...,*) 



determine the non-Gaussianity of the steady state distri- 
bution in position of the oscillator. As an example we 
give the results for some of them 

r (4) = 3Aq;4 

*> 4Af 4 w 2 7(3 7 2 + 4o; 2 ) 

C (5) = 24A« 5 



5M 5 w 2 (247 4 + 58 7 2 w 2 + 9w 4 ) 
5Aa 6 



9 " 3A/ 6 w 2 7 (107 4 + 37 7 2 u; 2 + 16w 4 ) 



(16) 



As already seen for the second and the third order cu- 
mulants above, the even order cumulants depend sensi- 
tively on the friction strength (oc I/7) while the odd cu- 
mulants saturate in the weak dissipation regime. Note, 
however, that friction must be always finite to reach the 
steady state, otherwise the oscillator will continuously 



heat up. Hence, for 7 <C ujq the even higher order cu- 
mulants always dominate against the odd ones, whereas 
in the overdamped limit 7 3> loq all cumulants are sup- 
pressed with Cq^ cx l/^ N ~ l \ As a consequence, the 
symmetry of the force field acts like a filter such as to 
favor even order cumulants with the tendency to es- 
tablish a symmetric non-Gaussian steady state distri- 
bution. For instance, the relative weight of the non- 
Gaussian cumulants with respect to the Gaussian one 

9n = [C^ /cfJ]{Mu) /a) n - 2 reads for weak friction: 
93 = (7/wo)/6,ff4 = 3/8,. 9 5 = ( 7 M))48/45, g 6 = 5/24. 



V. FREE PARTICLE WITH DRUDE DAMPING 

The above findings for a harmonic system can also be 
used to describe the dynamics of a free Brownian par- 
ticle subject to Drude damping and driven by Poisso- 
nian noise. In this case, the damping kernel is 7(4) = 
jlod exp(— uint) which in the limit of the Drudc frequency 
cj£> — > 00 reduces to the white noise situation. The most 
interesting regime is that of small lud though. Namely, 
the Langevin equation 



q{t) + / ds^t - s)q(s) = + r,(t)}/M 



(17) 



can be mapped to that of an harmonic oscillator by 
putting z(t) = J Q dsj(t — s)q(s). One finds 

z + u D z + 7 z = 7w B [£(i) + V(t)]/M (18) 

with q(t) = q(0) - f * ds [Mz(s) + £(s) + r]{s)]/M. Now, 
the Drude frequency takes the role of an effective friction 
constant, while 7 corresponds to the bare frequency of the 
fictitious z-oscillator. Accordingly, as in (|4| the inhomo- 
geneous part of the solution is determined by the response 
function ([3]) with 7 — > up and ui — > uj = \Ju)d1 ~ w £>/4- 
For free Brownian motion the most interesting quan- 
tity is the velocity q(t) and its cumulants which can now 
be inferred from the findings of the previous sections. 
The resulting expressions are rather complex, but ap- 
proach a stationary state after an oscillatory transient 
period (cf. Fig. [T]). While for larger ojd/j the third mo- 
ment (q{t) 3 ) quickly saturates, large amplitude oscilla- 
tions appear on a long-lasting transient period of time 
for small Drude frequencies. Correspondingly, the z- 
oscillator turns from an overdamped to an undcrdamped 
motion. In the stationary limit one finds 



lim (q{tf 



{2M 1 k B T + Xa 2 ) 1 



UJ D 



2jlo d 



lim (q(t) 3 ) = Xrt 



57 + 2u>d 
37 2 + 67W 



D 



(19) 



A realization of this model is given by a voltage-biased 
tunnel junction in series with an impedance consist- 
ing of an i?C-circuit (Ohmic resistance R + capacitor 
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C) and subject to Poissonian noise. In this case, the 
phase across the junction is related to the voltage via 
4>{t) = (e/h) f*dsV(s) and deviations = — eVt/H 
correspond to the generalized coordinate q in (fTTf Q. 
Its conjugate variable is the charge Q. The Drudc fre- 
quency is identical to the inverse of the i?C-time of the 
impedance. Hence, the velocity correlations (fl~9|) spec- 
ify charge fluctuations at the junction due to the non- 
Gaussian noise. According to Fig. [1] large charge oscilla- 
tions occur on long time scales for high impedances (large 
i?C-time). 
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FIG. 1. Dynamics of the third moment (q{t) 3 ) in units of 
Aq 3 /7 of a free particle immersed in a heat bath with Drude 
damping and subject to Poissonian noise in units of 7. The 
Drude frequency decreases from bottom to top: cjd/7 = 1 
(black, solid), u)d/~i = 0.1 (blue, dashed), ojd/j = 0.02 (red, 
dotted). 



VI. QUANTUM MECHANICS IN 
COORDINATE SPACE 



In the remainder we will generalize the results of 
Sec. HVI to the quantum regime, first in coordinate space 
and then in energy space representation. In coordinate 
space the most straightforward way is to follow the lines 
described in [ToL [Tl| as long as we arc mainly interested 
in steady state expressions. The linearity of Heisen- 
berg's equation of motion for the coordinate operator q(t) 
means that the quantum response function coincides with 
the classical one. Accordingly, the driving by Poissonian 
noise is taken into account exactly as in Qj. As a conse- 
quence, the non-Gaussian noise cumulants are identical 
to the results specified in (fT4")) and (fT!)|) . Only the sec- 
ond cumulant carries quantum mechanical information 
via the damping dependent equilibrium variance due to 



a heat bath with ohmic spectral density, i.e. 



(A 
?o 2 



7T(A H 
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9X + 

47T 
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6>A 

47T 

(20) 

where = ujqHP and X± = (7/wo) ± \J (7/wo) 2 ~ 4. Fur- 
ther the bare ground state variance is <?q = h/MQ and 
^ denotes the di-Gamma function. 

The reduced density operator of a harmonic system 
is obtained exactly within the path integral formalism 
by tracing out the Gaussian heat bath (for details see 
[111 ]) and then averaging over the noise drive, i.e., p(t) = 
(Trjj{VF(i)})^ where W(t) is the density operator of the 
full compound (system-l- reservoir) . In coordinate repre- 
sentation p(q,q') = (q\p\q') the steady state takes the 
form 



p st (x,r) = pm (r + -\p(t)\r - -) 



[T-q v {t- 



■(21) 



with r = (q + q')/2, x = (q — q') and q n as in &■ The 
momentum variance (p 2 }/3 reads as given in El, its 
specific form is not of relevance for the following discus- 
sion. The above distribution is no longer Gaussian with 
respect to the mean coordinate r but rather attains an 
asymmetry and long tails due to the Poissonian noise 
drive. However, it is still normalized according to the 
bare Gaussian normalization factor. For a weak noise 
amplitude a, an expansion of the diagonal part yields 



p st (x,r) = pfs{x,r) 



1 



(I 2 ) 



9»j 



3(9 2 }/ 



6{q% 



(22) 



(n) 

with the classical cumulants Cq and the bare thermal 
equilibrium pp(x, r) [results from (p?Tj) by putting q, q =0]. 
The asymmetry of the coordinate distribution /93(g) = 
[pst(0, q) - Ast(0, -q)]/pff(0, q) can now be probed to re- 
trieve in leading order the skewness (third order cumu- 
lant) of the noise drive. Its magnitude depends crucially 
on the dimensionless factor K3 = ((q 2 ) 3 p/qo)(^ + 2j 2 /u!q) 
in the nominator of ps such that at a fixed coordinate a 
small K3 yields a large P3. According to (fl?0"]) . for fixed 
damping K3 decreases with decreasing temperature. Par- 
ticularly intriguing is thus the limit of zero temperature, 
where one has 



1 2 o 



2^; 



[l-fArctan(^)] 7/^0 < 2 



T (A+-A-) ln ( A +/ A -) 7/^o > 2 



(23) 



with lu = \J uJq — 7 2 /4 as above. In the domain of weak 
friction, this leads to K3 « 1/8 — 37/(87rwo), while for very 
strong friction k 3 rj 16wo[ln(7/wo)] 3 /(7r 3 7) (cf. Fig. 
Hence, ps(q) is most pronounced for low damping at 
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7/wo ~ 0.25 where K3 attains a minimum, or in the over- 
damped limit where K3 decreases with increasing j/ujq. 

Recent experiments with current biased Josephson 
junctions as threshold detectors operate in the clas- 
sical regime. Then, for moderate to large friction the 
contribution of the third cumulant to the escape rate fol- 
lows up to a numerical factor from the asymmetry of the 
coordinate (superconducting phase) distribution in the 
harmonic well of a tilted washboard potential @- Our 
above results may provide signatures for the impact of 
quantum fluctuations when formally evaluating the con- 
tribution of the skewness at the position of the barrier 
top. Since for fixed friction, K3 decreases with decreasing 
temperature the relative impact of non-Gaussian noise is 
enhanced due to quantum effects. This strongly indicates 
the advantage for operating in the low temperature do- 
main which, however, necessitates a consistent quantum 
mechanical description of the full anharmonic problem. 




FIG. 2. Dimensionless factor K3 the size of which deter- 
mines the magnitude of the asymmetry in the coordinate 
distribution P3(q) vs. friction at zero temperature (see text 
for details). The inset displays a blow-up for weak damping 
with the dashed line corresponding to the undamped value 
^3(7/^0 = 0). 



bath H-q- An additional 'reservoir' £Tng is added which 
generates non-Gaussian noise. A general formulation is 
then based on a Hamiltonian 



H = H*+Hn + H 



I,B 



H 



I,NG 



H 



NG- 



(24) 



The interaction is bilinear in the reservoir and system 
degrees of freedom, i.e., 

■Hi, b + #1, ng = 9bQ ^2 c i x i + -9ng-R M , (25) 



where Q and R are system operators, the thermal bath 
Hb is modeled as a collection of harmonic oscillator de- 
grees of freedom Xj, and M denotes an operator of the 
non-Gaussian reservoir. We do not need to specify the 
non-Gaussian reservoir H^g (which in general may be 
very complicated) since only its correlation functions en- 
ter the formulation as we will sec below. The reser- 
voirs are statistically independent with CjXj)-Q = 

(M) NG = (MJ2j w)b,ng = 0. 

The goal is now to derive an approximate equation 
of motion for the reduced density p(t) within the Born- 
Markov perturbation theory up to third order in the cou- 
pling constants <?b and g^G (HI- After a straightfor- 
ward calculation and using the eigenstates of the system 
Hs\n) = E n \n), one obtains in the interaction picture 



(m\p s ,i(t) \n) 



kj 



(l\PsAt)\k) 



2 2 * 3 

9b dG 1 #NG dNG _ ^NG ^NG 
^2 knml ' knml ^3 ^knml 



(26) 



Here, the second order contributions contain conven- 
tional Redfield tensors in rotating wave approximation 
the explicit form of which is well-known in the literature 
fl2j . The non- Gaussian reservoir gives rise to a third or- 
der term C^^ ml that is lengthy and not very illuminating. 
Important to note is though that it carries system matrix 
elements of the form 



<n| {a±a?f \ri) 



(n| a±a^ \l) (l\ a±a^ \m) (m\ a±a^ \n') 



VII. QUANTUM MECHANICS IN ENERGY 
SPACE 



In the previous section we have demonstrated the pres- 
ence of non-Gaussian cumulants in the steady state distri- 
bution in coordinate space already in third order. Here, 
we focus on corresponding properties in energy space. 
These could in principle be inferred from the former re- 
sults, however, we want to put our analysis in a somewhat 
broader context and consider a full dynamical equation 
(master equation). 

For this purpose, one follows the conventional pertur- 
bativc approach for a system H$ weakly coupled to a heat 



(27) 

These matrix elements vanish if n = n' which has direct 
consequences on the appearance of this non-Gaussian 
contribution in the steady state. 

To see that in detail, we consider a two level sys- 
tem H§ = hwa z with the interaction Hamiltonian H\ = 
Cx(9ngM + 9B^2j c j x j)- The <r x ,a z denote Pauli ma- 
trices where |0), |1) is the eigenbasis of a z . At very low 
temperatures this system approximates the harmonic os- 
cillator sufficiently well. From ([21))) the relevant equations 
of motion in the Schrodinger picture are obtained as 



Poo 
Poi 
Poi 



-(WV 
2A/900 - 



h Wbi) Poo - ■+ 
A - 2£>'poi + (2D" 



W 10 
h u)p' 01 



(28) 
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where we exploited that pu + poo = 1 and poi = Poi 
ipQi = pIq- The conventional transition rates read 



IT" 



10/01 



ds 



&{e(s)£(0)) + ^(M(s)M(0)) 



(29) 



with £ = ^2jCjXj. The other second order contribution 
D = D' + iD" introduces a Lamb-shift and is given by 



D 



ds 



ik({£( s ),£(0)}} + ?ML({M(s),M(0)}) 



(30) 

with {, } being the anti-commutator. In (|2"5|) the leading 
non-Gaussian contribution appears as the real- valued pa- 
rameter 



according to (|26p off-diagonal elements p s t,n^m are zero 
and one can concentrate on the populations p s t,nn- The 
part in (f22j) originating from the third cumulant is anti- 
symmetric in r, while the bare equilibrium pp(x,r) is 
symmetric with respect to x,r —> —x,—r. The same 
is true for the product of the two wave functions with 
n = m so that the integral vanishes and the population 
is indeed independent of the third cumulant. 

We note that in contrast to the finding of the 
weak damping limit, the result (|22p is valid for any 
coupling to the thermal bath and only perturbative 
in the non-Gaussian noise. For stronger dissipation 
off-diagonal elements remain finite Pst.n^m 

^ and the 

third cumulant appears in the steady state distribution. 



A 



h 3 



ds / du(M(s)M(u)M(0)+M(0)M(s)M(u)) , 
Jo 

(31) 



It is derived from the tensor C^^ nl by assuming noise 
frequencies to be much larger than the system frequency 
uj (white noise limit). 

Given specific initial conditions, the above master 
equation can now easily be solved. For the ground state 
population the result is 



Poo(t) 



-rt 



r 

4Ap ' 1 (0) _ 



Poo(0) 



rt 



e 2Tt cos ojRt - 1 



(32) 



using the abbreviations T = W w + Wqi and = 
2uD" + uj 2 — T 2 /4. Obviously, the third moment in- 
fluences only the transient dynamics but is absent in the 
steady state behavior which here coincides with the ther- 
mal distribution. The reason for that can be traced back 
to the structure ([27)) : in steady state detailed balance for 
the populations dictates this non-Gaussian contribution 
to be absent. 

This can also be seen from the result in coordinate 
space (fT3|) directly. The energy representation follows 
from 

Pst.ram = (n\p st \m) 

dqdq'%l) n (r + p st (r, x) ip m (r - |^33) 

with wave functions i/) n (<l) = ( n \q) ancl eigenstates \n) of 
the harmonic system. Now, in the weak damping limit 



VIII. CONCLUSIONS 

In this paper we considered the paradigmatic model of 
a linear system interacting with a heat bath and subject 
to random noise with non-Gaussian characteristics. De- 
spite its simplicity, the steady state distribution of the 
relevant degree of freedom reveals a non-trivial depen- 
dence on parameters such as friction strength and tem- 
perature. In the classical regime this model also applies 
to charge fluctuations of a tunnel junction in series with 
a resistor+capacitor-impedance as heat bath. The re- 
duced density operator in the quantum case displays an 
enhanced sensitivity to non-Gaussian noise properties in 
coordinate representation at lower temperatures. Due 
to the symmetry of the force field in coordinate space, 
however, these properties are absent in steady state dis- 
tributions in energy space. Our findings are of relevance 
for all kinds of linear systems and may be used as a start- 
ing point for generalizations to nonlinear detectors such 
as Josephson junctions. 
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